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RESEARCH  SUMMARY 


Fall  frost  injury  to  159  seedlots  of  2-year-old 
Douglas-fir  seedlings  growing  at  a  nursery  in  Coeur 
d'Alene,  ID,  provided  an  opportunity  to  verify 
experimental  models  that  describe  genetic  variation  in 
freezing  tolerance  to  populations  from  the  Inland 
Northwest.  Experimental  results  of  three  previous 
studies  were  updated  and  then  were  combined  to  pro- 
vide a  single  model  of  freezing  tolerance.  According  to 
the  model,  elevation  and  geographic  location  of  the 
seed  source  accounted  for  67  percent  of  the  variance 
in  freezing  injury  to  167  populations.  Injury  predicted 
by  this  model  to  nursery-grown  seedlots  was  signifi- 
cantly correlated  (r  =  0.67)  with  the  injury  observed  in 
the  field. 
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INTRODUCTION 

In  the  Inland  Northwest,  Douglas-fir  (Pseudotsuga 
memiesii  var.  glauca)  occupies  a  broad  range  of  environ- 
ments that  v£iry  in  frost-free  period  from  60  to  150  days 
(USDA  Forest  Service  1965).  Populations,  however,  are 
physiologically  attuned  to  only  a  segment  of  this 
environmental  gradient  (Rehfeldt  1984).  When  acclimat- 
ing in  a  common  environment,  populations  within  north- 
ern Idaho  (Rehfeldt  1979a,  1979b),  western  Montana 
(Rehfeldt  1982),  and  central  Idaho  (Rehfeldt  1983)  dis- 
play a  differential  tolerance  to  early  faU  frosts.  Mathe- 
matical models  developed  for  each  of  these  three  geo- 
graphic regions  have  related  genetic  variation  to  the 
elevation  and  geographic  location  of  the  seed  source. 
These  models  have  generally  shown  that  populations 
from  areas  characterized  by  short  growing  seasons  suffer 
less  injury  from  fall  freezing  than  do  populations  from 
mild  environments. 

An  opportunity  to  verify  models  was  provided  by  fall 
frost  injury  to  2-year-old  Douglas-fir  seedlings  growing 
at  the  USDA  Forest  Service  nursery  at  Coeur  d'Alene, 
ID.  This  opportunity  also  provided  the  stimulus  for 
(1)  updating  eill  three  models  by  using  recent  procedures 
for  constructing  adaptive  landscapes  and  (2)  developing 
a  single  model  to  describe  genetic  variation  within  the 
InlEind  Northwest. 


MATERIALS  AND  METHODS 

Three  studies  (table  1)  of  tolerance  to  fall  freezing 
assessed  differentiation  of  Rocky  Mountain  populations 
in  common  gardens  but  under  much  different  experimen- 
tal conditions.  All  studies  involved  a  different  set  of 
populations  acclimating  in  different  years.  In  one  study 
tolerance  was  assessed  throughout  acclimation,  but  in 
the  other  two  tolerance  was  observed  on  a  single  date. 
In  one,  a  natural  frost  induced  differential  injury  to 
seedlings  in  a  field  test;  in  the  other  two,  a  laboratory 
chamber  was  used  to  freeze  detached  twigs.  Injury  from 
the  natural  freeze  occurred  on  a  single  night  and 
represents  a  single  treatment  that  induced  a  mean  injury 
of  only  8  percent;  that  in  the  laboratory  resulted  from 
numerous  treatments  that  produced  a  range  in  mean 
injury  from  nearly  zero  to  almost  100  percent.  As  a 
result,  the  freezing  tolerance  exhibited  by  a  population 
in  a  given  study  was  not  directly  comparable  to  the 
tolerance  of  populations  in  other  studies. 

Consequently,  the  development  of  a  model  describing 
genetic  variation  in  freezing  toleremce  for  the  entire 
region  required,  first,  a  single  data  set  normally  dis- 
tributed about  a  common  mean  and,  second,  a  data  set 
within  which  the  ranking  of  populations  accurately 
depicted  genetic  variation  uninfluenced  by  the  different 
testing  procedures.  While  the  former  condition  was 


Table  1— Summary  of  studies  of  freezing  tolerance  to  Douglas-fir  populations 


Characteristic 


Study  1^ 


Study  2^ 


Study  3^ 


Geographic  region 
Number  of  populations 
Origin  of  freeze 
Time  of  freeze 
Mean  injury  (pet) 
Binomial  variance 
Experimental  model 

Independent  variables 

R2 


Northern  Idaho 
43 

Laboratory 
Sept.  to  Nov. 
49 
0.0058 

7 
0.38 


Western  Montana 
56 

Laboratory 
mid-September 
39 
0.0042 

6 

0.35 


Central  Idaho 
73 
Nature 
mid-September 
8 

0.0010 

16 
0.48 


^Rehfeldt  1979b. 
^Rehfeldt  1982. 
^Rehfeldt  11983. 
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solved  with  transformations,  the  latter  was  met  by  scal- 
ing data  sets  according  to  the  differential  performance  of 
populations  common  among  studies. 

As  a  first  step,  the  effect  of  differential  mean  injury 
for  the  studies  was  eUminated  by  defining  a  deviation: 

where  Py  is  the  proportion  of  injured  seedlings  from 
population  i  in  study  j,  and     is  the  mean  for  study  j. 
The  constant  was  added  to  ensure  that  P '  >  0  and  to 
provide  deviates  about  a  mean  value  of  50  percent. 

Next,  normality  was  obtained  by  means  of  the  logistic 
transformation  sviited  for  binomial  data  (Finney  1971).  A 
logit,  Y;.,  was  weighted  according  to  the  variance  in  each 
study: 

Y;.  =  ln[(l  -  P.:)/P.:]/[(P.)(1  -  P.)nj-i] 

where  n.  is  the  number  of  populations  in  study  j. 

Finally,  all  data  were  scaled  to  those  from  northern 
Idaho  according  to  the  differential  performance  of  popu- 
lations common  among  studies.  Thus,  Y/was  defined 
such  that: 

1.  If  j  =  i,y:=y,, 

2.  If  j  =  2,  Y/=  Y.2  +  SF^ 

3.  If  j  =  3,  Y/  =  Yi3  +  SF^ 

where  the  scaling  factor,  SF^,  is  the  average  difference 
between  test  1  and  test  2  performance  of  two  popula- 
tions common  to  studies  1  and  2  and  SF^  is  the  average 
difference  between  test  1  and  test  3  performance  of  three 
populations  common  to  studies  1  and  3. 

These  procedures  produced  a  single  data  set  of  freez- 
ing injury  to  167  populations  that  was  normally  dis- 
tributed about  a  common  mean.  To  combine  data  sets  in 
this  manner,  however,  does  not  allow  the  effectiveness  of 
the  procedures  to  be  tested  statistically.  Effectiveness 
can  only  be  interpreted  according  to  the  degree  by  which 
the  combined  model  (1)  coordinates  the  individual  models 
and  (2)  is  biologically  plausible. 

The  combined  model  of  genetic  variation  in  freezing 
tolerance  was  developed  from  a  stepwise  regression  pro- 
gram for  maximizing      (SAS  1979).  The  program  used 
elevation,  latitude,  longitude,  northwest  departure,  and 
southwest  departure  plus  their  squares  as  independent 
variables.  The  northwest  and  southwest  departures  were 
obtained  by  rotating  the  grid  of  latitude  and  longitude 
by  45°.  Geographic  variables  were  nested  within  regions 
to  allow  a  given  independent  variable  an  opportunity  to 
describe  contrasting  patterns  of  variation  in  different 
geographic  regions.  Consequently,  26  independent  varia- 
bles were  available  for  stepwise  regression  of  the  general 
form: 

Y.'  =  /?„  +  )3,E.  Ef  +  E  7  X     +  X2 

where 

E;  is  the  elevation  of  population  i 

Xy^  is  geographic  variable  k  for  population  i  in  region  j 

/3q,  )3j,  iSg,  7  jk-  aiid      are  regression  coefficients 

k  =  1  .  .  .  4 

j  =  1  .  .  .  3. 

Adequacy  of  the  stepwise  model  was  determined  by 
the  level  of  significance,  the  size  of  the  residual  mean 


square,  and  a  lack  of  patterning  displayed  by  the 
residuals  (Draper  and  Smith  1981).  The  model  meeting 
these  criteria  was  used  to  predict  frost  injury  to  the 
seedlots  represented  in  nursery  production  beds. 

Frost  injury  to  2-year-old  seedlings  growing  at  the 
Coeur  d'Alene  Nursery  occurred  on  September  28,  1984, 
when  temperatures  reached  —  3  °C  after  a  prolonged 
period  of  warm  weather.  Differing  injury  occurred  to  179 
seedlots,  which  comprised  the  5  million  2-year-old 
Douglas-fir  seedlings  being  grown.  Of  these  seedlots,  159 
were  from  the  regions  represented  by  the  experimental 
model. 

In  late  October,  four  plots  of  100  seedlings  were  estab- 
lished at  equidistant  intervals  within  each  seedlot, 
regardless  of  size.  The  sample  thus  comprised  about  1.5 
percent  of  the  planted  trees.  Injury  was  scored  according 
to  browning  of  leaves  and  was  recorded  as  a  proportion 
of  injured  trees  for  each  seedlot.  After  the  logit  transfor- 
mation, these  binomial  data  were  then  correlated  with 
those  predicted  by  the  model  that  had  been  developed 
from  experimental  data. 

RESULTS  AND  DISCUSSION 

The  best  fitting  stepwise  regression  model  included  16 
variables  and  accounted  for  67  percent  of  the  variance  in 
freezing  tolerance  of  populations.  The  model  describes 
geographic  cUnes  for  a  common  elevation  (fig.  1)  and 
elevational  cUnes  for  various  geographic  localities  (fig.  2). 
Percentage  scales  in  these  figures  were  converted  from 
the  weighted  logit  transformation  by  using  a  variance 
among  populations  that  was  pooled  for  the  three  studies. 

Geographic  patterns  portrayed  by  the  model  for  a  con- 
stant elevation  indicate  that  populations  from  north- 
central  Idaho  exhibit  the  most  injury  (smallest  logits). 
In  fact,  mean  values  for  populations  from  this  area  devi- 
ate from  the  mean  by  more  than  two  standard  devia- 
tions (about  18  percent  injury).  From  this  area,  freezing 
toler^ce  of  populations  increases  in  all  directions,  but 
clines  are  particularly  steep  across  the  rugged  Bitterroot 
Range  to  the  west  and  across  the  Salmon  Drainage  to 
the  south. 

Geographic  cLines  at  a  constant  elevation  follow  geo- 
graphic patterns  of  frost-free  period  (U.S.  Department  of 
Commerce  1968)  even  though  geographic  trends  in  frost- 
free  period  are  not  independent  of  elevation.  Thus,  the 
populations  from  north-central  Idaho  that  exhibit 
highest  freezing  injury  (lowest  weighted  logit)  come  from 
an  area  tj^jically  having  long  frost-free  periods  (120  to 
150  days  at  800  m  elevation).  And,  corresponding  to  the 
genetic  patterns  of  figure  1,  the  average  frost-free  season 
gradually  decUnes  to  the  north  (90  days  at  900  m)  and 
rapidly  drops  to  less  than  60  days  in  the  rugged  moun- 
tains to  the  east  and  southeast. 

The  elevational  clines  (fig.  2)  are  presented  for  locali- 
ties on  three  of  the  geographic  isopleths  (fig.  1).  Like  the 
geographic  clines,  elevational  clines  (fig.  2)  closely  parallel 
the  declining  frost-free  period  associated  with  increasing 
elevation.  According  to  Baker  (1944),  the  frost-free 
period  changes  by  about  90  days  across  an  elevational 
interval  of  1,000  m.  In.  association,  regression  models 
predict  that  populations  from  800  and  1,800  m  at  the 
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Figure  1— Geographic  patterns  of  genetic  variation 
presented  as  isopleths  of  equal  performance  for  the  mean 
elevation  (1,420  m).  Isopleths  are  arranged  relative  to  the 
mean  predicted  injury  (x)  with  an  interval  betv^/een  isopleths 
equal  to  the  standard  deviation  which  corresponds  to 
approximately  9  percent  injury.  Shading  represents  the  dis- 
tribution of  Douglas-fir,  and  dots  locate  populations 
represented  in  the  experimental  data. 


 I  I  I  I  I  I  I  I  I  I  I 

800  1600  2400 

ELEVATION  OF  SEED  SOURCE,  M 

Figure  2— Elevational  dines  of  freezing 
injury  for  localities  along  three  geographic 
isopleths  (figure  1).  ▲  =  northern  Idaho, 
•  =  western  Montana,  ■  =  central  Idaho. 
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same  geographic  locality  will  differ  in  freezing  injury  by 
about  25  percent  when  mean  injury  is  50  percent.  The 
model  also  describes  a  nonlinear  cLine  that  was  undoubt- 
edly accentuated  by  a  lack  of  injury  to  most  central 
Idaho  populations  from  above  1,600  m.  The  nonlinear 
trend,  however,  was  also  apparent  in  studies  for  western 
Montana  (fig.  2)  and,  therefore,  is  not  completely 
attributable  to  the  low  mean  injury  associated  with 
study  3.  Apparently  elevations  above  1,800  m  are  of 
relatively  uniform  severity,  and  additional  adaptive 
differentiation  fails  to  occur  in  Douglas-fir. 

By  describing  genetic  variation  in  relation  to  environ- 
mental gradients,  the  model  could  be  used  to  predict 
relative  frost  injury  to  the  seedlots  growing  in  the 
production  nursery  at  Coeiu-  d'Alene.  The  relationship 
between  predicted  injury  and  observed  injury  was  not 
only  statistically  significant  (1  percent  level)  but  also 
accounted  for  45  percent  (r  =  0.67)  of  the  variance  in 
observed  injury  (fig.  3).  There  is  nO  doubt,  therefore, 


that  the  model  was  effective  in  predicting  actual  injury. 
Effectiveness,  moreover,  is  not  compromised  by  the 
scatter  about  regression  (fig.  3).  Relatively  large 
residuals  should  be  expected  when  the  injury  observed 
at  Coeur  d'Alene  involved  (1)  a  low  mean  injury 
(7  percent),  (2)  trees  spread  across  8.5  acres  of  produc- 
tion beds,  and  (3)  a  nonrephcated,  systematic  planting  in 
which  seedlots  from  the  same  Nationail  Forest  were  in 
adjacent  beds. 

Results  support  three  general  conclusions.  First,  the 
observed  injury  at  Coeur  d'Alene  vahdates  the  model  of 
genetic  variation  developed  from  experimental  data. 
Second,  artificial  freezing  of  detached  twigs  mirrors  the 
relative  amounts  of  damage  expected  for  intact  plants 
subjected  to  natural  freezing.  And  finally,  relatively 
steep  elevational  and  geographic  clines  describe  genetic 
variation  in  freezing  tolerance  of  Douglas-fir  populations 
in  the  Rocky  Mountains. 
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Figure  3— Relationship  between  observed 
frost  damage  to  seedlots  at  the  Coeur 
d'Alene  Nursery  and  injury  predicted  by 
regression  models. 
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tolerance.  According  to  the  model,  elevation  and  geographic  location  of  the 
seed  source  accounted  for  67  percent  of  the  variance  in  freezing  injury  to  167 
populations.  Injury  predicted  by  this  model  to  nursery-grown  seedlots  was 
significantly  correlated  (r  =  0.67)  with  the  injury  observed  in  the  field. 
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INTERMOUNTAIN  RESEARCH  STATION 


The  Intermountain  Research  Station  provides  scientific  knowl- 
edge and  technology  to  improve  management,  protection,  and  use 
of  the  forests  and  rangelands  of  the  Intermountain  West.  Research 
is  designed  to  meet  the  needs  of  National  Forest  managers, 
Federal  and  State  agencies,  industry,  academic  institutions,  public 
and  private  organizations,  and  individuals.  Results  of  research  are 
made  available  through  publications,  symposia,  workshops,  training 
sessions,  and  personal  contacts. 

The  Intermountain  Research  Station  territory  includes  Montana, 
Idaho,  Utah,  Nevada,  and  western  Wyoming.  Eighty-five  percent  of 
the  lands  in  the  Station  area,  about  231  million  acres,  are  classified 
as.  forest  or  rangeland.  They  include  grasslands,  deserts,  shrub- 
lands,  alpine  areas,  and  forests.  They  provide  fiber  for  forest  in- 
dustries, minerals  and  fossil  fuels  for  energy  and  industrial  develop- 
ment, water  for  domestic  and  industrial  consumption,  forage  for 
livestock  and  wildlife,  and  recreation  opportunities  for  millions  of 
visitors. 

Several  Station  units  conduct  research  in  additional  western 
States,  or  have  missions  that  are  national  or  international  in  scope. 
Station  laboratories  are  located  in: 

Boise,  Idaho 

Bozeman,  Montana  (in  cooperation  with  Montana  State 
University) 

Logan,  Utah  (in  cooperation  with  Utah  State  University) 

Missoula,  Montana  (in  cooperation  with  the  University  of 
Montana) 

Moscow,  Idaho  (in  cooperation  with  the  University  of  Idaho) 
Ogden,  Utah 

Provo,  Utah  (in  cooperation  with  Brigham  Young  University) 
Reno,  Nevada  (in  cooperation  with  the  University  of  Nevada) 


